Gravitational wave constraints on the shape of neutron stars 
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We show that there is a direct relation between upper limits on (or potential future measurements 
of) the m = 2 quadrupole moments of slowly rotating neutron stars and the I = m = 2 deformation 
of the star's surface, in full general relativity, to first order in the perturbation. This relation only 
depends on the star's structure through its mass and radius. All one has to assume about the star's 
constituents is that the stress-energy tensor at its surface is that of a perfect fluid, which will be true 
with good accuracy in almost all the situations of interest, and that the magnetic field configuration 
there is force-free, which is likely to be a good approximation. We then apply this relation to the 
stars which have direct LIGO/ Virgo bounds on their I — m — 2 quadrupole moment, below the 
spin-down limit, and compare with the expected surface deformations due to rotation. In particular, 
we find that LIGO observations have constrained the Crab pulsar's I = m = 2 surface deformation 
to be smaller than its I — 2, m — deformation due to rotation, for all the causal equations of state 
we consider, a statement that would not have been able to be made just using the upper bounds on 
the I = m = 2 deformation from electromagnetic observations. 

PACS numbers: 04.30.Tv, 04.40.Dg, 97.60.Jd 
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I. INTRODUCTION 

In the absence of significant internal stresses, objects 
bound by gravity are highly symmetric on large scales. 
Indeed, Lindblom and Masood-ul-Alam [HQ have shown 
that all static perfect fluid relativistic stars are spheri- 
cal (with reasonable assumptions about the equation of 
state). And even if one includes rotation, then equilib- 
rium states of fluid stars must have a high degree of sym- 
metry, as shown by Lindblom : They must be axisym- 
metric, if one takes dissipation into account, and must 
have reflection symmetry across the equator, at least in 
the Newtonian approximation (though this is conjectured 
to hold for relativistic stars, as well). 

However, situations in which the internal stresses are 
nonnegligible are not uncommon. Indeed, it is obvious 
that many of the smaller (though still gravitationally 
bound) objects in the solar system have large-scale in- 
ternal stresses that are comparable to their self-gravity. 
This is particularly true for asteroids and the smaller 
moons, which can be highly nonsphcrical (see, e.g., Q 
and Fig. 2.2 in Q), but is even true for the larger moons 
and the terrestrial planets, to a less immediately appar- 
ent degree. (For instance, see the model for the Earth's 
gravitational field given in Q.) What is perhaps surpris- 
ing is that such massive, strong-gravity objects as neu- 
tron stars can also have large enough internal stresses 
to produce a nonnegligible large-scale deformation, with 
reasonable theoretical input. In particular, as first appre- 
ciated by Ruderman 0, the star's solid crust can support 
a large-scale deformation, while an internal magnetic field 
will necessarily deform the star, as was first noted by 
Chandrasekhar and Fermi . 

In most cases, these deformations will be quite small 
(maximum fractional deformations of ~ 10~^, with con- 
siderably smaller deformations in most realistic cases). 
But there are certain scenarios in which the deforma- 



tion could be (relatively) substantial (> 10"^, possibly 
as much as tens of percent in extreme cases): If there is a 
solid exotic phase (e.g., the hadron-quark mixed phase) 
in the star's core, the entire star is solid (as is the case 
for crystalline superconducting quark stars), or there is a 
large internal magnetic field of > 10^® G. (See, e.g., (9l.[l0j 
for the solid case and (lT| - [l3t for the magnetic case.) 

The mechanisms for creating large deformations on 
nonaccrcting neutron stars arc not clear, though it is 
worth pointing out that neutron stars are expected to 
be born with some deformation from the violent super- 
nova explosion that forms them, and it is possible that 
this deformation may be frozen in to the solid parts, and 
only relax slowly over time, as discussed in |14| . See 
also Sec. 2.2 in [l5| for a discussion of scenarios involving 
accretion. However, deformed neutron stars are an inter- 
esting enough prospect for gravitational wave detection 
that it is worthwhile considering the more extreme cases, 
and seeing what bounds there are on large deformations. 
In most cases, the best bounds on the large-scale de- 
formations of neutron stars come from electromagnetic 
observations of their spin-down, which (as discussed in, 
e.g., [l^) constrain the star's I = m = 2 quadrupole 
moment, albeit with a factor of ^ 2 uncertainty due to 
the star's unknown moment of inertia. There are also 
tighter, though less direct, bounds that can be obtained 
from electromagnetic observations by folding in the star's 
braking index 17 1 . 



But the cleanest bounds come from gravitational wave 
observations, and a prominent accomplishment of the 
first generation of large-scale laser interferometer grav- 
itational wave detectors has been setting direct upper 
bounds on the to = 2 quadrupole moment of cer- 
tain known neutron stars, below the bounds that are 
given by electromagnetic observations of the spin-down, 
or indirect limits given by the star's age [ij, IT8l - [20j . 
See [H, [2l|, [2^ for recent reviews. It is customary to 
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quote these upper bounds in terms of a fiducial elliptic- 
ity for the star. This fiducial ellipticity gives a direct 
measure of the star's deformation in the uniform den- 
sity, Newtonian limit, but it is a priori unclear how to 
relate this to a measure of the deformation in more re- 
alistic cases. It is our purpose here to show how one 
can convert such fiducial cllipticities into a more physi- 
cally meaningful, fully relativistic measure of the shape 
of the star's surface. (While there is also the possibility 
of an m = 1 quadrupole deformation, as discussed by 
Jones [23j . we will not consider this situation here, since 
there are no final results for searches for such radiation — 
though see [l^j for preliminary results. Similarly, we do 
not consider radiation from higher multipole deforma- 
tions, which is suppressed by factors of the star's rota- 
tional speed over the speed of light compared to the I = 2 
radiation, and thus not searched for.) 

The relativistic shape measure we define is inspired by 
the measure of the horizon deformation of a tidally dis- 
torted black hole used by Taylor and Poisson [23| and 
involves the ratio of the star's longest polar circumfer- 
ence to its equatorial circumference. (Similar ratios have 
also been used in numerical studies of distorted black 
holes [26l-l28|.) The conversion from fiducial ellipticity to 
relativistic shape measure only depends upon the star's 
structure through its mass and radius (with reasonable 
assumptions about the magnetic field and shear stresses 
at the star's surface). This simple conversion comes from 
the simple relation between the metric on the star's per- 
turbed surface and the surface value of the time-time 
component of the metric perturbation (given by Damour 
and Nagar [l^), and the relationship between the surface 
value of this component and the amplitude of the star's 
quadrupole gravitational radiation, given in [13 (the re- 
lation with the quadrupole moment was first given by 
Hinderer [lOl). (As discussed in [l^, these relations as- 
sume slow rotation, using the results of Ipser [3l| and 
Thorne dH to show that it is appropriate to calculate 
the gravitational waves from the rotating star using the 
quadrupole moment calculated in the static limit. This 
is quite a good approximation for the stars we consider, 
which are rotating well below the Kepler limit — at most 
- 10% of it.) 

One can compute the deformation of the star due to 
rotation in the same manner, using Hartle's classic cal- 
culation [s^ of the metric of a slowly rotating relativistic 
star. While there is no reason to assume any sort of cor- 
relation between a star's I ^ m ~ 2 deformation and its 
rotational deformation,^ the rotational deformation gives 
a convenient scale against which to compare the bounds 



^ Indeed, millisecond pulsars, which will have the largest rota- 
tional deformations, for a fixed mass, are constrained by their 
small spin-downs to have some of the smallest I = m = 2 
quadrupoles — see, e.g.. Table 1 in [Tj; one can obtain a star's 
spin-down ellipticity by dividing the ellipticity bound given in 
the table by the value of hl^'^"/hf. 



on the I = m = 2 deformation. We then give these con- 
versions (from fiducial ellipticity and rotational frequency 
to surface deformation) for a variety of equations of state 
(EOSs), including the case of a strange quark star (for 
which one needs to modify the calculations slightly). We 
also compare the bounds on the I ~ m = 2 surface de- 
formation with the surface deformation due to rotation 
for the four stars with known spin periods for which the 
LIGO/Virgo bounds on the I ^ m ~ 2 deformation are 
near or below the (fiducial) spin-down limit. (We also 
note that models of strange quark stars with large radii 
can give significantly larger moments of inertia than the 
maximum of 3 x 10^^ g cm^ assumed in the LIGO/Virgo 
papers, which increases the uncertainty on the spin-down 
limit.) Additionally, we show that these results will re- 
ceive negligible corrections from non-perfect fluid contri- 
butions in the cases of interest (with the possible excep- 
tion of corrections to the force- free nature of the magnetic 
field configuration, which could be large enough to affect 
the results). 

The paper is structured as follows: We define the shape 
measure, give some intuition about its properties, and 
detail how to compute it from the fiducial ellipticity in 
Sec. ini We then compute the surface deformation due 
to rotation in Sec. IIIIl and compare the sizes of these 
two deformations for a variety of equations of state and 
pulsars of interest in Sec. IIVI Finally, we conclude and 
consider the outlook for future bounds on the I = m ~ 2 
deformations of known neutron stars in Sec. |Vl In the 
Appendix, we show that the corrections to these results 
due to the non-perfect fluid nature of real neutron stars 
are likely negligible in realistic cases. 

II. THE RELATIVISTIC SHAPE MEASURE 
A. Definition 

First recall that the fiducial ellipticity is defined by 



where I^x and lyy are the star's actual moments of in- 
ertia about the axes other than its rotation axis (the 
difference between the two is related to the star's m = 2 
quadrupole moment, which is the observable quantity) 
and = lO'*^ g cm^ is the fiducial moment of inertia of 
the star about its rotational axis. [As is discussed in, e.g.. 
Sec. VI B of [l^ the actual moment of inertia of neutron 
stars is expected to lie in the range (1-3) x 10"*^ g cm^, 
on theoretical grounds, but there are no measurements 
of this quantity for any neutron star. And as we shall 
see, one can obtain even larger moments of inertia, up to 
at least ~ 5.5 x 10**^ g cm^, for models of strange quark 
stars with large radii.] 

In contrast, our relativistic measure of the deforma- 
tion will be constructed solely from quantities on the 
star's surface, which could be measured (in principle) 
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by a physical observer, viz., the star's equatorial cir- 
cumference, Soq, 22 1 and its longest polar circumference. 



Spol, 22, 



givmg 



esurf, 22 



Spol, 22, max ~ ^oq, 22 
Scq, 22 



(2) 



(cf. the calculation of the surface deformation of a tidally 
deformed black hole in Sec. VIII of Taylor and Pois- 
son [2^, and the calculations of the shape of deformed 
black holes in [26l - [28l |).^ One can compute the quanti- 
ties entering this definition using the metric on the star's 
surface, so all wc need to do is write this metric in terms 
of the star's quadrupolc moment (along with its unper- 
turbed mass and radius). 

To gain some intuition about the properties of this 
deformation measure, we consider a Newtonian ellipsoid, 
with semi- axes (1 + 'i])R, (1 — 'i])R, and R, for which 
we have esurf, 22 = '?/2 (to first order in 77). [Here we 
have used the expression 27ra[l — e^/4 + 0{e'^)] for the 
circumference of an ellipse with semi-major and semi- 
minor axes a and b and eccentricity e = ^yl — b^/a?.] 
If the ellipsoid has uniform density, then I^x ~ lyy = 
4r]I,, = (8/5)77Afi?2 [note that we are using the true 
moment of inertia about the z-axis, I^z = {2/5)MR^, 
not the fiducial one here, and still working to first order 
in 77], so its Newtonian ellipticity is eN,22 = 477, giving 

Esurf, 22 = ^N, 22/8. 



B. Computation 

Let us now relate the relativistic shape measure to 
the star's fiducial ellipticity. The star's surface metric 
is given in standard spherical coordinates by Eqs. (85) 
and (86) in Damour and Nagar [2^ 



ds^ = i?2 ( 1 + 2^r,„J [de"^ + sin^ Odcf)'^), (3) 



where R is the star's undeformed radius, Yi,„ is a spher- 
ical harmonic (we shall only consider Z = rri = 2 in our 
discussion), and the deformation is given by [Eq. (92) 



^ While we do not need this generality here, since we only con- 
sider I = 2, m a {0, 2} first-order deformations of a spherical 
star, note that one can define this surface ellipticity for an ar- 
bitrary two-dimensional manifold that is homeomorphic to the 
2-sphere, imbedded in M? , and possesses a preferred axis (here 
given by the axis of rotation) . The generalization of the lines of 
longitude are the geodesies passing through both intersections of 
the preferred axis with the manifold (i.e., the "poles"), and the 
equator can be defined as the locus of points which have equal 
geodesic distance from both poles. Also note that in the rota- 
tional case, these circumferences can no longer be determined 
directly from timelike or null geodesies, due to frame dragging. 
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Here C := 2GM/Rc? is the star's compactness, where M 
is the star's mass. (Note that Damour and Nagar define 
a compactness of half this, which they denote by c.) This 
expression assumes a nonrotating star, but this is quite a 
good approximation for the stars in question, since they 
are rotating slowly enough that one can also treat their 
rotation as a perturbation of a static star, and the two 
perturbations are independent to first order (since their 
multipolar structure is different ). The corrections due 
to rotation will thus be at most ~ 1%. Additionally, as 
discussed in [loj . the gravitational quadrupolc radiation 
from a slowly rotating star can be computed from the 
static quadrupolc moment of the star with no rotation, 
as shown by Ipser [3l| and Thorne [3^ . 

We now want to relate this to the I = m ~ 2 
quadrupolc moment of the deformed star. Here we con- 
sider the quadrupolc moment amplitude Q22 given by 
Eq. (1) in [l^l in the Newtonian limit. We thus sub- 
stitute the expression for Hq(R) in terms of Q22 from 
Eqs. (39) and (41) in [l^, giving 
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(5) 



/22(C), 



where we have defined Ho{R) to stand for Ho{R) 
with ci ^ 1, and also defined the "correction factor" 
/22(C). [Since Ho{R) only depends on C, the correc- 
tion factor indeed only depends on C. Also note that 
dlogi/o(i?)/rflogi? = '~d\ogHoiR)/d\ogC.] 

It is necessary to make a small change to this compu- 
tation to treat the case of strange quark stars, where the 
density does not go to zero at the surface. This affects 
the present calculation because the pressure perturbation 
no longer goes to zero at the surface, so one must add on 
— {G/c^)ttHo{R)R'^P^ to the expression for SR/R given 
in Eq. where p_ is the density just inside the star's 
surface. [This comes from inserting the expression for K 
from Eg. (39a) in Ipser [U into Eq. (88) in Damour and 
Nagar [23, using Eq. (35) in to relate the pressure 
perturbation to the density, and noting that Ipser's pi is 
the same as —6p, in the notation of [l^l- We neglect the 
shear stress terms in all of these relations, which we shall 
show is a good approximation in the Appendix.] How- 
ever, this only increases the star's deformation by < 5% 
in the case we consider. 

We now compute the desired arclengths, finding (to 
first order in 6R/R) 
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and 

Spol, 22 = 2i?o 
= 2ttRo 



15 SR. 
32^ ~R 



sin 6*008 200 



d0 



1 + 



15 SR 
1287r~R 



cos 2(/ic 



(7) 



where we have used Re Y22 = ^/T5j32n sin^ 6* cos 2(/) and 
considered the longitude line with azimuthal angle 
We can now write Csm-f, 22 in terms of Q22 , or the fiducial 
I = m = 2 ellipticity efid, 22 = \/^t^7^Q 22/ 1 zt ^ using 



Eqs. @ and ©, yielding 



/15^ c4 Q22, 15 c4 /fd 



512 G2 M3 • 

/22(M,C)efid, 22- 



64 G2 M3 • 



(8) 



Here f22{M,C) gives the conversion factor between the 
{I = m = 2) fiducial ellipticity and the I = m = 2 surface 
deformation. [Nota bene: Since /22(C) > 0, we have 
SR/ R < 0, and thus the maximum value of Spoi, 22 is 
given by taking cos20o = 

Note that in the Newtonian limit, we have esurf, 22 = 
{3/8){I^^/MR'^)efid, 22- This does not agree with the 
Newtonian calculation for the uniform density ellipsoid 
given at the end of Sec. Ill Al which gives a coefficient of 
5/16, but that is to be expected, since the present com- 
putation assumes a fluid star, whose surface will deform 
to be an cquipotcntial of its perturbed gravitational field, 
while the ellipsoid was taken to be rigid. 



III. CALCULATION OF THE ROTATIONAL 
DEFORMATION 

It is interesting to compare upper bounds on an ^ = 
m = 2 deformation to the ^ = 2, m = deformation 
that will be induced in the star due to its rotation. Here, 
one can use the slow- rotation results of Hartle [s^ to per- 
form the calculation. We start by recalling that [Eq. (86) 
in Damour and Nagar [2^; cf. Eq. (25a) in Hartle and 
Thorne [H 1 



SR Sr 1 , , 



(9) 



where Sr/r gives the fractional position of star's deformed 
surface, and K{R) is the (/ = 2, m = piece of the) an- 
gular component of the metric perturbation, evaluated 
at the star's surface."^ [Here, following Hartle, we take 
the angular dependence to be just the Legcndrc polyno- 
mial portion of Y20, viz., (Scos^fl — l)/2, without the 



^ We neglect the I = change to the star's radius, since it gives 
a second-order correction to the calculation of the relativistic 
shape measure. 



normalization factor of y^b/Air present in the spherical 
harmonic. Additionally, note that we have reversed the 
sign of Damour and Nagar 's K to correspond with the 
conventions of Hartle |33| and [3.] Eq. (146) in Har- 
tle [3^ gives an expression for Hartle's ^2/0, which is the 
same as our Sr/r, so we have 
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h2iR) = A 
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C2/2 



1 -C 



(11) 



[Eqs. (137) and (139) in Hartle 



, ri is the magni- 
tude of the star's angular velocity, and J = Izz^ is the 
magnitude of its angular momentum, where Izz is the 
star's (true) moment of inertia. Here we have noted that 
Hartle's C := 2/C - 1. 

The amplitude A is given by solving Eqs. (125) 
and (126) in Hartle, using the procedure given above 
Hartle's Eq. (146). However, note that in our experi- 
ence it is possible to streamline this procedure slightly: 
We integrate the inhomogeneous equations starting from 
values of at tq, the inner radius where we impose our 
inner boundary condition when solving the enthalpy ver- 
sion of the Oppenheimer-Volkov e quat ions j35| — see the 
discussion at the end of Sec. HI in — and take A ~ 1 
(this A is different from the amplitude) in the boundary 
conditions when solving the homogeneous equations. (A 
different value for A is equivalent to an overall scaling, 
by homogeneity.) We also convert the equations to en- 
thalpy form, as was first done for the frame-drag equation 
in Sec. 4.1 of K. Lockitch's thesis [l^. (The transformed 
frame-drag equation also appears in a slightly different 
form in Appendix A of [s^l-) One can then perform the 
matching at the surface using /i2 and its first derivative, 
and use the V2 matching as a check. (This check indicates 
that our calculations are accurate to better than ~ 2%, 
usually much better.) 

We also note that the surface deformation scales as fi^, 
so we can simply calculate Izz and A for a fixed central 
value of Hartle's frame-drag parameter w, and then scale 
our solutions to give the desired angular velo city , using 
the surface value of w and Eq. (47) in Hartle |33l. S imi- 
larly, from Eqs. (137) and (139)-(141) in Hartle [Hi, we 
have, noting that Hartle's v := ^2 + fc2 , and that K = 2k2 
[cf. Hartle's Eq. (66) and Eq. (24) in [13 ], 



A 



8G2M4 



C 



3(| 



1 log(l-C) 



(12) 



(1+C) 



Finally, we need consider the changes that need to be 
made to this calculation to treat strange quark stars, with 
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their nonzero surface density. Due to Hartle's method 
of locating the surface of the rotating star using a first 
integral to the equations of hydrostatic equilibrium, no 
change is necessary in that portion of the calculation. 
However, we do need to modify the surface matching 
used to obtain the amplitude. The calculation of the 
moment of inertia remains unchanged, since only the 
second derivative of w is discontinuous at the surface. 
But the contribution to Hartle's (j^)_R (the subscript R 
denotes a radial derivative, in Hartle's notation) from 
the density evaluated just inside the star's surface is 
-87r(G/c2)i?p /(I ~ C) [from Eqs. (26), (28), (29), and 
(40) in Hartle [33| ]. Thus, the matching of the solutions 
to Hartle's Eqs. (125) and (126) at the surface needs to 
be adjusted using the replacements 



3 c4 



^2, out ^ ^2, out 



1-Cj 1-C 
47ri?2p„/i2(i?) 47r G /2 C 



(13a) 



M 



R^p^uj^iR) 
1~C ' 



3 c4 V C 1-C 



(13b) 



where the subscript "out" denotes the solution outside 
the star and primes denote radial derivatives. 

One can compare the values for the quadrupole ob- 
tained using the Hartle slow-rotation expansion with the 
fits to a fully relativistic calculation for a 1.4 M© star 
given in Eq. (73) and Table 5 in Frieben and Rezzolla [l3| . 
Here one uses Eq. (26) in [s^l to obtain the quadrupole 
from the slow-rotation calculation (the expression given 
in [33| is in error). We find quite close agreement for 
the four EOSs we both consider — at worst, ^ 2% for the 
BBB2 and GNH3 EOSs, and at best, better than 0.1%, 
for the APR and SLy EOSs. (What we caU the SLy 
EOS is the same as Frieben and RezzoUa's SLy4; we give 
further discussion of these EOSs in Sec. |lVl) Our calcu- 
lations of the moment of inertia all agree to better than 
0.1%. 

Now doing the calculations of the equatorial and polar 
circumferences, we have, as previously, 
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(14) 



Eqs. ©-(US]). Thus, we have 

Scq, 20 ~ ^poL 20 _ 3 / dR 



Csurt, 20 



■Scq, 20 



4 V i? 



(16) 



to first order in 6R/R, where we have used the known 
scaling of the surface deformation with ft^ in the final 
equality to write Cgurf, 20 in terms of the star's spin fre- 
quency, i^, defining the conversion factor g2o{M,C). We 
have also introduced a minus sign, compared with the 
1 — 171 — 2 version, so that esurf, 20 will be positive. 

For the purposes of comparison, we relate our esurf, 20 
to the coordinate radius elligticity or flattening [used by, 
e.g., Frieben and Rezzolla [ij] — see their Eq. (68)] in the 
Euclidean limit (and to first order in the deformation). 
This coordinate radius ellipticity is defined by 



(17) 



where re and are the star's equatorial and polar radii, 
respectively, so we have Ec = 2esurf. 20 in the given limit. 
(Here we have noted that the eccentricity of an ellipse is 
related to the flattening by = ec/2 



el.) 



IV. RESULTS AND DISCUSSION 

Here we show the results of our calculations for a rep- 
resentative selection of EOSs, whose mass-radius curves 
and moments of inertia versus mass are illustrated in 
Fig. [TJ (All the calculations presented here were carried 
out using Mathematica 7's default methods.) Specif- 
ically, we have chosen the four EOSs from the lorene 
library [11] that are compatible with the Demorest et 
al. [3^ measurement of a 1.97 ± 0.04 Mp neutron star 
(the nucleonic EOSs APR g^, BBB2 HJ, and SLy [13 
and the hyperonic EOS GNH3 [H; note that BBB2 is 
only compatible within 2a). We also consider the nu- 
cleonic WFFl EOS [11], the most compact EOS consid- 
ered in (33, obtained from the RNS website [45|.* For 
more exotic EOSs, we use three of the hadron-quark 
hybrid EOSs from [i^ (Hyl, LKRl, and generic) and 
the strange quark matter EOS constructed in [l^] using 
the results of Kurkela, Romatschke, and Vuorinen |47| 
(KRVl). 



Spol, 20 — 2i?o / 
JO 



2nRo 



1 + 1(^4) (3cos20-i; 
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20 



1 /6R 



4 \ R 



20. 



de 



(15) 



where {dR/R)2o is calculated by combining together 



* Note, however, that the high compactness of stars obtained with 
the WFFl EOS is in part attributable to the fact that this EOS 
becomes acausal (i.e., has a sound speed greater than the speed of 
light) at densities well below the central density of the maximum 
mass stable stars. The same is true for the EOS we consider that 
generates the second most compact stars, the APR EOS, which is 
computed using a variational method, like the WFFl EOS. Thus, 
in our plots, we will mark the point at which stars constructed 
with these EOSs start to contain acausal matter. 
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FIG. 1: The mass-radius relation (left) and moment of inertia versus mass (right) for the EOSs considered in this paper. For 
the APR and WFF EOSs, the asterisks mark the maximum mass for which the stars do not contain any acausal matter (i.e., 
for which the central density is smaller than the density at which the EOS has a sound speed greater than the speed of light). 



We show the conversion factor f22{M,C) for these 
EOSs in Fig. [2] Note also that if one considers a canon- 
ical M = 1.4Afo, R — 10 km neutron star, then its com- 
pactness is - 0.42, and /22(1.4M0, 0.42) ~ 0.22. Thus, 
we see that the fiducial ellipticity generally gives a rea- 
sonable impression of the size of the star's surface de- 
formation. Of course, the actual size of the deformation 
depends upon the star's undeformed size: One can give a 
dimensionful measure of the surface deformation by mul- 
tiplying Esurf by the star's radius, R (recalling that the 
star's equatorial circumference is unchanged by the de- 
formation, to linear order). We show this for the LIGO 
upper limits on the Crab's quadrupole moment in Fig. [5J 
ta king Cfid, 22 = 10~^ (around the best upper limit given 
in |14{). And while we have used the current upper limit 
on the Crab's fiducial ellipticity, one can scale these re- 
sults linearly to apply to any other fiducial ellipticity (for 
the Crab or any other neutron star). 

Note, however, that if one considers elastic deforma- 
tions, then the SLy EOS would not be able to support 
a quadrupole large enough to create a fiducial elliptic- 
ity anywhere close to lO"**, unless one assumes a very 
low-mass star — see Figs. 5 and 6 in and Fig. 3 in 
Horowitz . This is simply because the crust is the 
only solid portion of neutron stars with nucleonic EOSs, 
like the SLy EOS (and even hyperonic EOSs), and the 
crustal shear modulus is (relatively) small, compared 
with the shear moduli of the higher-density solids that 
may be present in more exotic EOSs. Thus, while there 
do not yet exist calculations for the maximum deforma- 
tions that could be sustained by the other hadronic or 
hyperonic EOSs, they also should not be able to support 
such large quadrupoles, with standard crustal composi- 
tions, since stars constructed with these EOSs have com- 
pactnesses that are roughly the same as those obtained 
using the SLy EOS. And even the three hybrid EOSs 
shown would only be able to support a quadrupole large 
enough to create such a large fiducial ellipticity for high 



masses — see Fig. 9 in [lO|. However, as shown in Fig. 11 
in [l3|, one could easily obtain such a large quadrupole 
with the KRVl EOS, assuming that the strange matter 
is in a crystalline state, with a shear modulus given by 
the Mannarelli, Rajagopal, and Sharma [i^ calculation. 
And strong enough internal magnetic fields could sustain 
such a large quadrupole for all EOSs. See, in particular, 
the calculations of magnetic deformations by Frieben and 
Rezzolla [l^ for four of the nucleonic EOSs wc consider. 

Now considering the rotational deformations, we show 
the conversion factor g2o{M, C) and the dimensionful sur- 
face deformation for the Crab in Fig. [2] Again, one can 
scale the results for the Crab to any other (slowly rotat- 
ing) pulsar, using the scaling of the surface deformation 
with v^. 

The behavior of all these quantities is as expected: As 
the star becomes more massive and more compact, one 
finds that a given I = m = 2 quadrupole moment (i.e., 
a given fiducial I = m = 2 ellipticity) translates into 
a smaller surface deformation, as one would anticipate, 
since the star's moment of inertia is increasing. This 
trend reverses for the highest masses, however, where 
the relativistic suppressions of the quadrupole due to 
the increasing compactness (discussed in Sec. V of [l0| ) 
now dominate. Additionally, the more massive, compact 
stars do not deform as easily, and this is seen in the 
mass dependence of the rotational deformation. In the 
strange quark star case, one even sees that the decrease 
in deformability with increasing mass is strong enough to 
dominate the increase in radius. In general, for a fixed 
rotation rate and fiducial ellipticity, the I = m = 2 sur- 
face deformation is the largest and the rotational defor- 
mation is the smallest for the most compact stars. This 
is illustrated using the bounds and expectations for the 
dimensionful surface deformation of the Crab pulsar and 
the most compact EOSs we consider (all nucleonic) in 
Fig. El 

In Table m wc compare the LIGO/ Virgo upper bounds 
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FIG. 2: The I = m = 2 (upper left) and rotational (lower left) conversion functions [/22(M, C) and g2o{M,C)\ and the maximum 
I — m — 2 (upper right) and expected rotational (lower right) surface deformations for the Crab pulsar in dimensionful form 
(i.e., Resuri, 22 and Resuii, 20), all versus mass and for the various EOSs we consider. For the maximum I = m = 2 Crab pulsar 
surface deformations, we use the LIGO upper bound on the fiducial ellipticity of £^"^^22 = 

10"*. Also, for the APR and WFFl 

EOSs, the asterisks mark the maximum masses that do not contain acausal matter. 



on the I = m = 2 deformation with the expected rota- 
tional deformations for the four neutron stars for which 
the LIGO/ Virgo bounds are lower than indirect limits, 
and there is a measured spin period (viz., the Crab and 
Vela pulsars, J0537-6910, and J1952+3252) [H,!!!. We 
do this assuming the fiducial 1.4 neutron star and 
considering all the EOSs used in this paper. But recall 
that there is no reason to assume any correlation be- 
tween the I = m = 2 and I ~ 2, m ~ deformations. We 
merely quote the numbers for rotation as a convenient 
scale against which to compare the upper bounds on the 
nonaxisymmctric deformations. 

Additionally, note that the last two stars we consider 
require a moment of inertia larger than the fiducial one 
(of 10"*^ g cm^), or a significantly closer distance than the 
estimate used in [3l in order for the LIGO upper bound 
to beat the spin-down limit. This is not too much of a 
concern for J0537-6910, as one only needs a factor of ^ 2 
larger moment of inertia, which is easily provided by any 
number of EOSs (see, e.g., this paper's Fig. [T] and Fig. 6 
in [m). However, for J1952-I-3252, Q claims that LIGO 
likely just misses beating the spin-down limit, even given 



the uncertainties in distance and moment of inertia. We 
have still included this star here, since the requisite mo- 
ment of inertia of > 4x 10*^ g cm^ would be easily attain- 
able with the KRVl EOS for higher masses.^ Addition- 
ally, there is a fair amount of uncertainty in this pulsar's 
distance: Using H i measurements, Verbiest et al. (5^ 
give distance bounds of 3±2 kpc, so on the small side, this 
is significantly closer than the distance of 2.5 kpc used 
in [l4| . and would be enough to put the LIGO bound 
under the spin-down limit, even assuming the fiducial 
moment of inertia. (However, note that reducing the 
distance would reduce the inferred limit on the fiducial 



In considering the effects of the moment of inertia on the spin- 
down limit, be aware of the following typographical errors in Il4l : 
In that paper's Eqs. (1) and (7), /sg should be ^33'^ and 



^38 ' 



respectively. The first of these errors is also present in an inline 
equation in the first paragraph of Sec. 1 of [iSl and Eq. (14) 
of I20II [though Eq. (1.2) in II6II is correct], and the second is 
also present in the first paragraph of Sec. 3 of [T^ . though the 
analogous Eqs. (6.1) in II6II and (15) in |20|| are correct. 
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TABLE I: The maximum I — m = 2 and expected rotational (Z = 2, m = 0) surface deformations (both dimensionless and 
dimensionful) for the four pulsars for which there are LIGO/ Virgo upper bounds below the spin-down limit, for the various 
EOSs we consider, assuming a 1.4 Mq star. We also give the stars' rotational frequencies i/, and the upper bounds on fiducial 
ellipticity from the LIGO/Virgo results, e^^^22- The values for are valid for the years 2005-2007, corresponding to the time 
of the gravitational wave observations. (While the Vela pulsar observations were later, the Vela pulsar's spin-down is relatively 
slow, so its spin frequency is not affected, to the given accuracy.) Due to the Crab pulsar's large spin-down, its current spin 
frequency is 29.7 Hz [SO], which reduces its rotational surface deformation by ~ 0.7% compared with the values given in the 
table. All values for the deformations are rounded to two significant figures. 
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FIG. 3: The maximum I = m = 2 (lower curves) and ex- 
pected rotational {I = 2, m = 0; upper curves) dimensionful 
deformations (i.e., i?esurt, 22 and Rtsuri, 20) of the Crab pul- 
sar versus mass for the hadronic EOSs we consider, using the 
LIGO bound of 

^^^^22 = 10 As before, asterisks denote 
the maximum masses that do not contain acausal matter for 
the APR and WFFl EOSs. 



ellipticity from the LIGO results, in addition to increas- 
ing the spin-down limit, while increasing the moment of 
inertia does not affect the limit on the fiducial elliptic- 



ity.) But note that even these large strange quark star 
moments of inertia are not quite large enough to make 
the LIGO hmit lie below a reasonable spin-down limit for 
J1913-I-1011, the pulsar with the next closest LIGO limit 
to its fiducial spin-down limit, without assuming that 
the pulsar is significantly closer than thought: The LIGO 
Umit for J1913-I-1011 is 4.9 times its spin-down limit, and 
the spin-down limit scales as lH'^ — cf. Eq. (1.2) in [l6| . 

Gonsidcring the EOSs presented here, we can see 
that gravitational wave observations have constrained the 
large-scale (i.e., / — 2) surface deformation on the Crab 
to be due to rotation, not to any sort of Z = m = 2 distor- 
tion. [The only exception comes from the near-maximum 
mass stars given by the WFFl EOS, with its high maxi- 
mum compactness of 0.668; note that such stars contain 
a significant region (~ 15% by volume) of acausal mat- 
ter.] This conclusion could not have been made solely 
from electromagnetic observations, since the star's spin- 
down only puts a limit on the I = m = 2 deformation 
that is at least ~ 7 times greater than the LIGO limit 
(and possibly a factor of a few more, given the uncertain- 
ties in the moment of inertia) — see Table 3 in [TJ]- This 
would only be sufficient to constrain the Crab pulsar's 
Z = m = 2 surface deformation to be smaller than its 
rotational surface deformation for stars with large radii. 
(Note that some interpretations of X-ray observations of 
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neutron stars (53|, imply that their radii are ~ 11- 
12 km, so the EOSs that give stars with larger radii 
may be disfavored, though these interpretations are by 
no means conclusive.) And even the smaller, less direct 
upper limit obtained by Palomba [Ol by folding in the 
star's braking index only bounds the deformation to be 
^ 3 times larger than the LIGO bound, which would be 
insufficient to ensure that the I = m ^ 2 surface defor- 
mation is smaller than the rotational surface deformation 
for massive, compact stars, as shown in Fig. [31 

Since the Vela pulsar spins at less than half the fre- 
quency of the Crab pulsar, and has an upper bound on its 
fiducial ellipticity that is an order of magnitude greater, 
the observations do not bound its / = to = 2 surface 
deformation to be smaller than its rotational surface de- 
formation, even for the case in which the two deforma- 
tions are the most similar, the LKRl EOS and a mass 
of ^ 1.8 Mq. Conversely, since J0537-6910 has a rota- 
tional frequency that is ^ 2 times the Crab's and the 
same bound on its fiducial ellipticity, we know that its 
rotational deformation is larger than its ? = to = 2 de- 
formation for all masses and EOSs, even for the most 
massive, compact stars with the WFFl EOS. (Indeed, 
the spin is large enough that this conclusion could have 
been obtained from the electromagnetic observations of 
the spin-down, even with the uncertainties in the moment 
of inertia.) The bounds on J1952-t-3252's I ^ m = 2 
surface deformation compared with its rotational defor- 
mation are similar to the Crab's, though somewhat less 
good, due to J1952-l-3252's somewhat smaller spin fre- 
quency and the larger upper bound on its fiducial elliptic- 
ity. In particular, the bound on J1952-|-3252's I = m ~ 2 
surface deformation will be larger than its rotational sur- 
face deformation for masses slightly larger than 1.4 
for the WFFl EOS, and for aU the other nuclconic EOSs 
for larger masses. 

Finally, for a close-to-home scale against which to com- 
pare these surface deformations, the Earth's (dimension- 
less) rotational surface deformation is ~ 1.7 x lO^'^, and 
the (dimensionless) I = m = 2 surface deformation of the 
geoid (i.e., of a gravitational equipotcntial near sea level) 
is ^ 2.8 X 10~^. These values were obtained using the 
value of the flattening of the Earth's reference ellipsoid 
and the I ^ m = 2 component of the Earth's gravita- 
tional field, both from the WGS84 version of the Earth 
Gravitational Model EGM2008 Q, along with the rela- 
tion between the flattening and the rotational surface de- 
formation given at the end of Sec. IIIIl The Crab pulsar's 
and J1952+3252's (dimensionless) rotational surface de- 
formations arc thus very similar to the Earth's for the 
EOSs with larger radii, if the pulsars' masses are close to 
1.4 M0. However, all of the bounds on the I = m = 2 
(dimensionless) surface deformation for the pulsars we 
consider are larger than the Earth's I ~ m = 2 (dimen- 
sionless) surface deformation, though the bounds (elec- 
tromagnetic or gravitational wave) for many of the faster 
rotating pulsars arc well below this (sec Table 1 in [l3|). 



V. CONCLUSIONS AND OUTLOOK 

We have shown that one can convert the bounds on 
the fiducial ellipticity of neutron stars (from gravitational 
wave or electromagnetic observations) into bounds on the 
shape of the star's surface, in full general relativity, to 
first order in the perturbation (and in the slow-rotation 
limit). We have given the conversion for a variety of 
EOSs, some purely nucleonic and others containing ex- 
otica, including the strange quark star case (with the 
requisite changes to the calculation) , and compared with 
the shape of the star's surface due to rotation. Here we 
find that the gravitational wave observations have con- 
strained the Crab pulsar's I = m = 2 surface deformation 
to be smaller than than its Z = 2, m = surface deforma- 
tion due to rotation, for all the causal EOSs we consider. 

For the other three pulsars for which LIGO/ Virgo ob- 
servations have beaten the spin-down limit, within the 
uncertainties on the star's distance and moment of iner- 
tia, we find that the Vela pulsar's relatively slow rotation 
means that the bound on its Z = to = 2 surface defor- 
mation is well above its rotational surface deformation, 
for all the EOSs we consider. On the other hand, for 
the more quickly rotating J0537-6910, the I = m = 2 
surface deformation is known to be below the rotational 
deformation from the electromagnetic observations of its 
spin-down. (Again, this holds for all the EOSs we con- 
sider.) For J1952+3252, the bound on the Z = to = 2 
surface deformation is below the rotational surface de- 
formation for all cases except for heavy stars constructed 
with the hadronic EOSs. 

Looking into the future, the prospects for improving 
these bounds are good. Indeed, as mentioned in Sec. 6 
of [lOl , Virgo-|- science run 4 (VSR4) was expected to be 
sensitive to fiducial cUipticitics of a fcwx 10~^ for the Vela 
pulsar and ^ 10~^ for the Crab pulsar and J1952+3252 
if it achieved its sensitivity goal, and preliminary results 
using data from this run indeed give improved bounds 
for the Crab and Vela pulsars [2J|. These improvements 
would bound the I = m ~ 2 surface deformations for the 
Vela pulsar to be below its rotational surface deforma- 
tion for EOSs with larger radii, and place the bounds for 
the Crab pulsar and J1952-I-3252 well below their rota- 
tional surface deformation for all masses and EOSs. (Of 
course, in the most optimistic scenario, one would obtain 
a measurement of the I = m = 2 surface deformation 
instead of a bound.) There are also several other pulsars 
for which the Enhanced LIGO and Virgo-f observations 
already made could be able to beat the spindown limit 
(as discussed in Sec. 6 of p^). 

Advanced LIGO and Virgo are expected to have sensi- 
tivities more than 10 times better than the initial inter- 
ferometers for the pulsars considered in , which would 
allow them to further improve the bounds on the pulsars 
already considered and to place bounds below the spin- 
down limit for many more pulsars (at least 47, based on 
the results in Sec. 2.2 of Pitkin |2l|). And looking even 
further into the future, Pitkin [21| predicts that the Ein- 
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stein Telescope will plaee bounds below the spin-down 
limit for hundreds of pulsars. Additionally, as shown in 
Fig. 3 of Pitkin [2lj, the majority of the stars for which 
these detectors are expected to beat the spin-down limit 
have (relatively) small spin frequencies, so constraining 
the I = m = 2 surface deformation to be smaller than 
the expected I — 2^ m = rotational deformation gives 
another (vaguely physical) benchmark for searches, be- 
low the spin-down limit. However, we must stress once 
again that there is no expected relation between the two 
deformations. 

On the theoretical side, it would be interesting to com- 
pute the rotational surface deformation of more rapidly 
rotating stars using either lorene (SS,] or rns ^] and 
compare with the slow-rotation predictions given here. 
Such a calculation, or the more involved calculations 
likely needed to obtain the bounds on the / = m = 2 
surface deformations of more rapidly rotating stars from 
bounds on their gravitational wave emission, could be of 
interest in interpreting results from next-generation de- 
tectors, which are expected to be able to beat the spin- 
down limit (by a bit) for one or two pulsars with spin 
frequencies of ^ 500 Hz (see Fig. 3 in Pitkin ^l|). How- 
ever, these frequencies are still low enough (compared 
to the Kepler frequency) that our slow-rotation approxi- 
mation is likely still reasonably accurate. (See, e.g., the 
small deviations from linear dependence of the spin pa- 
rameter on frequency through half the Kepler frequency 
in Figs. 4 and 7 in Lo and Lin [55|, and the values for the 
Kepler frequency in that paper's Fig. 2.) 

Likely more important would be the inclusion of the 
magnetic field in the calculation of the conversion from 
the I = m = 2 quadrupole moment to the surface de- 
formation: As discussed in the Appendix, it is possible 
that small departures from a force-free configuration at 
the star's surface could affect the conversion for stan- 
dard pulsar magnetic fields. And even if the magnetic 
field configuration is force-free to better accuracy than 
we need, the order-of-magnitude estimates of its effects 
on the computation given in Eq. (jA.ip suggest that for 
magnetar-level surface magnetic fields of ~ 10^^ G, and 
the fiducial ellipticities of around 10~^ associated with 
internal fields of this magnitude (as given in, e.g.. Sec. 7 
of Frieben and Rezzolla [l^), the effects of the magnetic 
field could be large enough to affect the results we have 
given for the shape. However, it is worth noting that the 
internal field could be significantly larger than the sur- 
face field (as discussed in, e.g., Corsi and Owen [11]), in 
which case the corrections would be significantly smaller. 
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Appendix: Corrections due to non-perfect fluid 
contributions 



Since we have obtained our results for the I = m = 2 
shape deformation using expressions that were derived 
assuming an unmagnetized perfect fluid star, we need to 
check that they are still valid in the cases in which we are 
interested, where the star is perturbed not by the tidal 
field of a companion (as in Damour and Nagar [1^), but 
by some internal stresses. (Since we assume that we can 
treat these internal stresses as a first-order perturbation, 
we do not need to worry about their effects on the ro- 
tational deformation, since the two perturbations will be 
independent, to the order we are computing.) It turns 
out that we can indeed apply their expression with no 
changes, if one assumes that the stresses at the star's 
surface (times GB? /c^ or 1/ p-c^) are much smaller than 
Hq{R), which we shall sec is the case in the majority 
of situations of interest. Explicitly, one needs to make 
this assumption in obtaining the expression for the star's 
perturbed surface using its enthalpy perturbation [see 
Eq. (26) in Damour and Nagar [2^ ], and also in writing 
the metric function K{R) in terms of Ho{R) and Hq{R). 

In the first case, the expression for the enthalpy used 
to obtain the position of the perturbed surface is ob- 
tained from stress-energy conservation. The relevant 
equation including the stress terms is given in Eq. (A2) of 
Ipser 



3l| . and with slightly different notation in Eq. (35) 



of |10|. In the second case, Ipser gives the relation be- 
tween K and Hq including the stress terms in Eq. (39a). 

We now consider the extent to which realistic sur- 
face stresses affect these results. Looking at the second 
case first, the star's magnetic field produces a stress of 
~ B^/Stt, so the fractional corrections from these stresses 
to the relation between K and Hq will be on the order of 



GB'^R^ 



G^B^M^R'^ 



c^Ha{R) nc^Q22Ho{R) 

3 / n \ 2 

10-9 



M 



1.4 M, 







R 



10 km 



B 



f 10-"* 



1012 g; Vefid,22. 

(A.l) 



[We have omitted the dependence on the star's com- 
pactness from Ho{R) in the final expressions here and 
later, though it will further suppress these corrections 
for higher-mass stars.] 

Considering the corrections to the stress-energy conser- 
vation equation used to locate the star's surface, we have 
fractional corrections on the order of RJ°'Fab/c^ PsHq{R), 
since the addition to stress-energy conservation from the 
electromagnetic field is V'^T^^ = J''Fab, where J° and 
Fab denote the (4-)currcnt density and Faraday tensor, 
respectively. [Compare Eq. (35) in (T^]; the magnetic 
field term would appear on the right-hand side of that 
equation, and the factor of R comes from the factor of 
VaYim that multiplies everything else in that expression.] 

We expect J''Fab (which gives the Lorentz force den- 
sity) to be small, compared to Ps/R, since the magnetic 
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field in the magnetosphere is expected to be close to force- 
free (as is discussed, e.g., at the end of Sec. 2.2 in [st}). 
If one takes the magnetosphere to be exactly force-free 
(as is frequently done in models of magnetically deformed 
neutron stars), then J°'Fab is zero. However, in the situ- 
ations we consider, J°'Fab could be of the same order as 
the Hq term, since in this case, the Lorentz force could 
help balance the gravitational perturbation, instead the 
pressure perturbation alone balancing it, which is what 
occurs in the absence of additional stresses. But treat- 
ing this properly is a task for far more detailed neutron 
star modeling (analogous to the issue with surface cur- 
rents discussed in Sec. II C of Corsi and Owen (s^), so we 
simply take the magnetic field configuration to be exactly 
force- free here. 

For the Crab pulsar, with a surface field of ^ 4x 10^^ G, 
one would have to consider a fiducial ellipticity of at most 
^ 10~^^ for the first of these corrections to approach the 
percent level. This is well below the predicted sensitiv- 
ity of even the proposed Einstein Telescope (see, e.g.. 
Fig. 3 in Pitkin 12111 ). and around the minimum ellip- 
ticity quoted in [18[ as being produced by an internal 
field of the same magnitude as the Crab's external field. 
(The Frieben and RezzoUa [l^ fits predict an ellipticity 
of ^ 10^^ for an internal field of that strength.) The 
other three pulsars considered in Ta ble II have external 
magnetic fields of similar magnitudes [58| , so this analy- 
sis holds for them, as well. 

On the other hand, if we consider a solid strange quark 
star, the surface shear stresses could be quite large, if the 
gap parameter does not decrease too much as one ap- 
proaches the star's surface and the high breaking strain 
assumed for the high-pressure interior continues to the 
surface — see the discussion in Sec. IV C of [l3l. How- 
ever, it seems likely that the stresses will be considerably 
smaller in any sort of reasonable case, since one expects 
both the gap parameter and breaking strain to decrease 
as the pressure decreases. And if one assumes that the 
surface is not strained much more than the star's average 
strain, then the surface stress cannot be too large in the 
cases of interest. 

Specifically, considering the KRVl strange quark EOS, 



we have ratios of (in order of magnitude, and putting 
in the ^ 15.5 km radius appropriate for a 1.4 Mq star 
with this EOS, in addition to the associated values of the 
Mannarelli, Rajagopal, and Sharma [i^ shear modulus) 

~^ V 1.41/0/ V15.5 kmj V 1033 erg cm-3 J 

"(—)■ 

V^fid, 22/ 

(A.2) 

and 

/x_cr_ _ G^M^H-a^ 
^ Ho{R)p-c^ ^ nc^Q22Ho{R)p- 

\iamqJ \ 10-2 ; Vefid,22y ' 

where p-d- is the shear stress just inside the star's sur- 
face (written as a product of shear modulus fj.- and shear 
strain <J-). (These come, as before, from the corrections 
to the relation between K and Hq and to the enthalpy 
expression used to locate the star's surface.) 

While it is possible that (T_ could be as large as ^ 10"^, 
around the breaking strain for the neutron star crust ob- 
tained in [59I |63| (and thus orders of magnitude larger 
than Efid, 22 for the stars we are considering), this high 
breaking strain comes from high pressure, so the break- 
ing strain at the surface, where the pressure goes to zero, 
will likely be much less. Moreover, if one assumes that 
the star's surface is strained at about the same level as 
the star's average strain, then we will have (t_ ~ efid. 22 
(cf. Fig. 11 in jlOf . remembering that one can linearly 
scale the maximum quadrupoles given there for a uni- 
form strain of 10"^ to any uniform strain). Thus, we 
feel comfortable quoting values for the strange quark star 
case, since the caveats seem fairly mild in reasonable sit- 
uations. 
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